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SUMMARY 


Measurements of actual runway roughness obtained by a pro file -survey 
method (engineer’s level) are presented. Data were obtained from a survey 
of a relatively rough runway and a smooth runway. The results of this 
study are presented as roughness profiles of the runways surveyed and in 
the form of power spectra. 


INTRODUCTION 


The frequency of occurrence of large load applications in routine 
ground airplane operations has caused a growing concern in regard to the 
roughness of landing and taxiing surfaces. In order to obtain informa- 
tion on this problem, it was thought desirable to make detailed measure- 
ments of the roughness characteristics of actual runway surfaces . As 
one of the initial steps in this study, measurements were made of two 
runways available to the National Advisory Committee for Aeronautics at 
Langley Field, Va. The two runways selected were known to be of very 
different degrees of roughness; one runway was considered relatively 
smooth whereas the other was considered rather rough - possibly rough 
enough to preclude active use . The measurements made are presented here 
directly as elevation profiles. In addition, the power spectra of the 
runway elevations were also determined and are presented in order to 
permit a description of the frequency characteristics of the runway 
roughness . 


SYMBOLS 


D distance over which moving average is taken 

L wavelength, ft 

m number of uniformly spaced points over the frequency range at 

which power estimates are desired 
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n number of equally spaced, elevations taken over the runway 

x distance, ft 

X arbitrary value of x, ft 

Ax space interval, ft 

y(x) random function of distance (runway height) 

R(X) autocorrelation function 

®(n) 


x 

n 

a 


SURVEY OF VARIATIONS IN RUNWAY HEIGHT 


A diagram of the landing and parking strips that presently exist 
at Langley Field is presented in figure 1 and shows the extent of the 
two surveys made . Runway 17-35 was chosen because it is considered 
representative of a satisfactorily smooth runway; the other runway chosen, 
12-30, is considered rough and is used only for parking* Both runways 
are of standard concrete construction. 

The roughness measurements were made by means of a surveyor's level, 
rod, and tape . This means was selected because it could be applied 
directly without the delays attendant on the development of special instru- 
mentation. In using this technique it was necessary to select an interval 
at which elevations would be obtained. It was thought that the frequency 
range between 0-5 and 35 cycles per second would be the region of prin- 
cipal concern for most airplanes . At a landing speed of 100 miles per 
hour this would correspond to wavelengths between about 300 and k feet. 

As a consequence, a reading interval of 2 feet was selected. This choice 
was dictated by two considerations: first, it was expected that there 

would be little variation in runway height at wavelengths less than k feet 
and, second, the communications sampling theorem (ref. 1), which states 
that sampling a disturbance at intervals of one -half the shortest wave- 
length present completely specifies the disturbance. 


power -spectral density function, defined by equation (l), 
ft 2 

radian/ft 

displacement distance, ft 

reduced frequency, 2 tt/l, radians/ft 


root-mean-square value of y; V/y^x) 
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Only 1,400 feet at the south end of runway 35 was surveyed so as 
not to interfere with ground traffic from runway 7- Visual observation 
suggested that except for the initial 100 feet, which is a sloping macadam 
overrun, this 1,400 feet is fairly representative of the remainder of the 
runway. Three thousand feet or nearly all of runway 12 was covered. The 
two surveys were conducted with greater ease than was originally expected. 
Approximately 6 hours was spent on runway 35 and slightly more than twice 
that time was spent on runway 12 . This means that on an average about 
115 readings per hour were obtained. The pace was steady but not hurried. 
Detailed runway elevations are given in table I, and figure 2 shows the 
runway profiles . These elevations are plotted about a zero arithmetic 
mean . 


POWER SPECTRA OF RUNWAY HEIGHT 
Definition of the Power Spectrum 


In addition to the actual runway -height profiles, it appeared desir- 
able, because of the random character of the height fluctuations, to deter- 
mine the power spectra of runway height. These power spectra would pro- 
vide a description of the frequency content of the runway height variations 
and be directly applicable to the calculation of airplane responses in the 
frequency plane . 

Since the runway roughness under consideration here is a space 
disturbance rather than a disturbance in time, it is desirable to define 
the power-spectral density function in terms of the frequency argument ft 
in radians per foot rather than the conventional argument c 0 in radians 
per second. In terms o f this frequency argument . the power -spectral 
density function of the disturbance y(x) is defined in the following 
manner ; 



1 

2jtX 



(1) 


where the bars indicate the modulus of the complex quantity. Equation (l) 
may be used to evaluate the power-spectral density function from observed 
data, but, in practice, the power -spectral density function may be deter- 
mined more conveniently and less tediously through use of a related func- 
tion, the autocorrelation function R(X), defined by 


R(X) 



( 2 ) 
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The autocorrelation function has the symmetrical property R(X) = R(-X) 
and is reciprocally related to the power-spectral density function by 
the Fourier cosine transformation in the following manner: 


®(«) 



R(X)cos ftX dX 


poo 

R(x) = / a>(n)cos nx cw 

J o 


From the foregoing it is seen that 


( 3 ) 


R(0) = y2(x) 


a 2 = f 4>(n)dfl 
^0 


( 4 ) 


where a is the root-mean-square value (or standard deviation) of the 
disturbance and is a convenient measure for a comparison of the overall 
roughness of the two runways surveyed. 


Evaluation of Power Spectra 

The actual evaluations of the power spectra were made by means of 
the numerical -calculation procedure derived by Tukey in reference 2 . 

This procedure is also described and discussed in reference 3 . 

As a preliminary to these power-spectrum calculations, it appeared 
desirable to make some modifications to the actual measured profiles . 
Examination of these profiles (fig. 2) indicated that the runway height 
exhibited large changes in elevation at very long wavelengths . These 
large changes at low frequencies have a tendency to complicate and distort 
power estimates at the higher frequencies because of the effective filter 
characteristics of the numerical estimators. Since, in the present study, 
there was little interest in the longer wavelengths, it was decided to 
avoid these adverse affects by removing the longer wavelengths . One 
simple and convenient way of prefiltering a disturbance is by the use 
of moving averages . This means was used in the present study to filter 
out some of the longer wavelengths . A moving average of the runway height 
for a 300-foot distance, as defined by 
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1 k=75 

y m ( x ) = ttz *>1 y(x + k Ax) 
k=-75 


was determined for each of the runways and is shown as the dashed curve 
in figure 2. The variations of runway height about this moving aver ag e 
were then determined and are shown in figure 3" As can be seen roughly 
from a comparison of figures 2 and 3, the main effect of this operation 
is the removal of the low-frequency components of the height variations . 

The actual filtering effects of this operation on the spectrum 
estimates are derived in detail in the appendix. The attenuation func- 
tion for the spectrum introduced by the moving average is given by 

(‘ ' f j 

where D is the distance over which the average is taken. The attenua- 
tion function is shown as a function of QD/2 in figure 4. A separate 
scale for wavelength L is also shown for the present value of 
D = 300 feet . It can be seen from the figure that the principal effect 
of the moving average in the present case is to attenuate the effects of 
the longer wavelengths, with the attenuation factor decreasing from 1 
at 300 feet to 0.490 at 400 feet and 0.132 at 600 feet. 

The actual steps involved in the numerical estimation of the spectra 
(ref. 2) are as follows: 

1. The autocorrelation coefficients were determined from the succes- 
sive values of elevation y lf y 2 , . . . y n according to the following 

numerical form of equation (2): 


R r 



n-r 


q+r 

q=l 


(r = 0, 1, 2, . . . m) (5) 


where = R(x) j X=r/«. 

2 . Initial or raw estimates of the power density were then deter- 
mined by use of the following numerical form of equation (3) : 
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Lr = 



m-1 

2 Rq cos 

q=l 


3^1 + IL. cos r«| 
m 111 * 


(6) 


where L_ = L(ft) : Q = — „ . These estimates have an effective filter 

1 m Ax 

which has the undesirable character of appreciable side-band areas and 
thus permits a wide diffusion of power . 

3 . Final or "smoothed" estimates of power density, which are estimates 
based on a more desirable and sharper filter, were then determined from 


t>Q = 0.54Lq + 0.461)! 

= 0 »25Lj»_i + 0 + 0 •23f , 2'_j.! ^ 


G m = 0.461^.! + 0.541m 


(7) 


The values of <t r obtained in the foregoing manner can roughly be con- 
sidered to be estimates of the average power over the frequency interval 

.(. r . < n < t—iliL. The distances of 1,400 and 3,000 feet covered 

m Ax m Ax 

in the survey and the interval Ax of 2 feet lead, to values of n of 
701 and 1,501, respectively, for the two runways. Evaluation was made 
with m = 40 . The autocorrelation coefficients obtained are shown in 
figure 5 and the resulting power estimates are plotted as a function of 
the reduced frequency 0 in figure 6 . Each point in the present case 
represents the average power in a frequency interval ± ^ about the 

value plotted. For clarity in presentation, some of the power estimates 
at the higher frequency are not shown but were used in obtaining the 
faired curves . A scale of wavelength L is also shown in figure 6 . 

The root mean square of runway height a, which is a convenient measure 
of the average roughness power, is listed in figure 6 for each of the 
runways in the figure . 

It may be recalled that the spectra of figure 6 represent the spectra 
of the runway height variations about the moving average . The effects 
of the moving -average operation on the derived spectra have been shown 
in the appendix to be equivalent to multiplying the actual spectra by the 
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attenuation function 


1 - 



In principle , then, it would appear 


that the effects of the moving -average operation could be removed by 
dividing the spectra of figure 6 by this attenuation function. However, 
in practice this cannot be done explicitly, since the spectrum estimates 
are averages of the power over finite band widths . A good approximation 
of the effect may be obtained, however, by dividing the faired spectrum 
estimates of figure 6 by the attenuation function and then averaging the 
resultant spectra over each of the band widths . It may be seen that these 
operations would have a negligible effect for the range of frequencies 
covered in f Igure 6 . 


DISCUSSION 


Examination of figures 2 and 3 indicates that the runway heights 
fluctuate in a random manner, with variations of as much as several inches 
for each runway. The profile for runway 12 on both figures shows wider 
fluctuations, particularly at the longer wavelengths. The overall varia- 
tions in height at the wavelengths of principal interest are best compared 
in figure 3; where the effects of the very long wavelengths (greater 
than 300 feet) have been largely removed by the moving average . It is 
clear from this figure that runway 12 is appreciably rougher than runway 35 > 
as was to be expected. The height profiles of figures 2 and 3j representing 
what may be considered as a satisfactory and an unsatisfactory runway, 
provide representative runway Inputs for response calculations . 

A more detailed comparison of the characteristics of the runway height 
variations is possible from the spectra of figure 6. The overall height 
fluctuations as represented by the root -me an -square values of 0 .057 feet 
for runway 12 and 0 .021 feet for runway 35 indicate that runway 12 is 
almost three times as rough as runway 35* By comparing the spectra it 
can be seen that runway 12 has 10 times the power of runway 35 at the 
longer wavelengths (300 to 500 feet) and about twice the power at the 
shorter wavelengths (below 50 feet) . The rapid decrease of power with 
increasing frequency displayed by both spectra is perhaps generally 
typical of runway height spectra and provides a guide to a representa- 
tive spectrum shape . Inasmuch as the spectra of figure 6 represent 
satisfactory and unsatisfactory levels of runway roughness, the heights 
of the spectra provide an initial guide toward the establishment of 
criteria for runway roughness . 
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CONCLUDING REMARKS 


As an initial step in the study of runway roughness, data were 
obtained from a survey of a relatively smooth runway and a runway which 
is considered rough. A surveyor's level, rod, and tape were used to 
obtain these data. In general, this method was found to be quick and 
relatively inexpensive . The results obtained are presented as profiles 
of runway height . The power spectra of runway height were also deter- 
mined and are presented. These results provide an initial guide toward 
establishment of criteria for runway roughnesses and axe suitable for 
airplane response calculations . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 1J, 195^ • 
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APPENDIX 

EFFECT OF MOVING AVERAGES ON THE POWER SPECTRUM OF A RANDOM FUNCTION 

Consider a random function f(x) . The centered moving average of 
f(x) over the interval D is defined by 


x + g 


S(X) = D / B f(x l )dx l 


(Al) 


The deviation of the disturbance about its moving average is then given 


D 
x+ — 


AAt — 

h(x) = f(x) - g(x) = f(x) - i J D 2 f( Xl )dX! 


(A2) 


'The spectral properties of h(x) as compared with those of f(x) are 
of concern; specifically, the effects on the spectrum of f(x) intro- 
duced by the operation of equation (A2) . 

The "finite" Fourier transform of equation (A2) may be vritten 


H(n) = / h(x)e -1 ^ x dx = 

J -X 


v -lux, 
f(x)e dx 


i r x r * 1 

D jvtj 


2 f(Ti)e” i ^ X d T l dx 


•X x- 


(A3) 

If the order of integration of the double integral term is inter- 
changed, there results 

H(Q) = J f(x)e -iax dx - i J D 2 J ^ f(Tj)e _inx dTi dx (Ak) 
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The x-integral of the second term may now be evaluated to give 


H(n) 



sin 



OD 

2 / 


|F(fl) 

(A5) 


where F(n) is the "finite" Fourier transform of f(x) 


F(Q) = f f(x)e '^ lX dx (A6) 

J -X 

The power spectrum of h(x) follows directly from equation (A5) in 
accordance with the following equation: 


® h (n) - lim ttW l R ( a )l 2 “ ^?H(Q)H(-Q) 

X >00 X >00 


(A7) 


Substitution of equation 
simple result: 


*h(n) = 


sin 

HD 

2 


(A5) into this equation yields the following 


lim -i- F(o)F(-Q) = 
X — 2nX 




$£ (fi) 


(A8) 


The squared term on the right side of this equation is the effective 
spectrum attenuation function or transfer function introduced by the 
moving-average process defined by equation (A2) . It is shown plotted as 
a function of f ID/2 in figure 4. Also shown in the figure is a scale 
of wavelength for D 3 300 feet, the value used in the present study. 

As a matter of added interest, an alternate derivation of equation (a 8) 
may be given which is slightly more involved but which is more general 
in the sense that the function g(x) does not have to be functionally 
related to f(x) . Equation ( Al) may also be written 


g(x) 


1 

D 



f(x 1 )k(x 


xi)dx x 


(A9) 
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where 


k(x) = 1 (-D/2 ^ x ^ D/2) 

k(x) = 0 (elsewhere) 


Now, from equation (A2) the power spectra of h(x) and g(x) can be 
shown to be related in the following manner (ref. 4) : 


®h(n) = ®f(fl) + ® g (a) - «f g (ft) - ® g f(fl) 


(A10) 


where <fp(fi) and $g(fl) are the power spectra of f(x) and g(x), 
respectively, and <ffg(ft) and $gf(fl) axe the cross spectra of f(x) 
and g(x) defined in the following manner: 


$>f g (a) 


® g f(n) 


lim F(-a)G(n) 
X — > 00 2nX 


> 


lim F(fl)G(-Q) 
X — > co 2nX 


(All) 


where F(u) is defined by equation (A6) and 


G(Q) = f g(x)e ^ X dx (A12) 

J -X 

The cross spectra in general have both real and imaginary terms, in con- 
trast to the simple power spectra which are always real. 

In order to determine the spectrum of h(x), it is necessary to 
evaluate each of the spectra on the right side of equation (AlO) . Con- 
sidering <J^;(n) first, it can be seen from equation (A9) that g(x) is 

related to f(x) by the conventional convolution or Duhamel integral 
and is thus obtained by a linear operation on f(x) . Using the simple 
relation between spectra of an input disturbance and an output response 
for linear systems, there is immediately obtained the following relation 
between g(x) and f(x) : 
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•g(0> 



(A13) 


where the squared, term represents the amplitude of the frequency- 
response function for the operation defined by equation (A9) . 

The cross-spectra terms may be evaluated in the following manner . 
From equations (A9) , (All), and (A12), for example. 


F(-ft) 


tf’fg(n) = lira 


r r - 

-X 


f(x^)k(x - xq)dxq^|e ^ X dx 
2*X 


(A14) 


In the limit, this equation reduces to 


(A15) 


The cross spectrum in this case has no imaginary term reflecting the 
lack of any introduction of a phase shift by the linear operation defined 
by equation (A2) . In the present case 

$fg(n) = (ai6) 

Substituting the results of equations (A13) , (A15) , and (Al6) into 
equation (A10) yields the same result as given by equation (A8) . 


® fg (n) = V n ) 


sin 


OD 

2 
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Distance along runway, ft 

Figure 2.- Variation of runway elevation about an a 
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(a) Runway 35* 



Distance along runway, ft 

Figure 3- _ Variation of runvay elevation about a 300~foot moving mean 










